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MathArt Fractal MUD 
Stewart Dickson (see [Ref1]) 


Stewart Dickson announces the unveiling of his lat- 
est interactive art work on the Internet: the MathArt 
Fractal MUD (Multi-User Dungeon). The piece is 
the multi-user exploration of a three-dimensional 
Julia set, which can be accessed using standard 
World-Wide Web browsers like NCSA Mosaic and 
Mosaic Communications’ Netscape. You can find it 
at the internet site whose URL (Uniform Resource 
Locator) is given in reference [Ref2] (see Refer- 
ences, below). 


Fractal MUD advances the theory that there is 
something metaphysical at the ultimate boundary 
of something we cannot fully comprehend. 
http://mathart.com also contains documenta- 
tion on creating physical sculpture of deterministic 
fractals. (See accompanying article, “Direct Sculp- 
ture of Deterministic Fractals”, below) 


The Julia set calculated in quaternion four-space is 
a graph of infinite complexity and detail. Human- 
kind never saw such a beast before 1982. We can- 
not view all four dimensions at once and we can 
never reach the finest level of detail in the graph. 
Throughout human history, the boundary of 
knowledge has been filled with superstition and 
spirituality. Therefore, the artist has placed in every 
chamber of Fractal MUD an Oracle which speaks 
from “Beyond the boundary”. 


The system uses the HyperText Markup Language 
(HTML) of the World-Wide Web as a program data 
stack to navigate the three-dimensional shadow of 
complex four-space. The user. sees a two-dimen- 
sional slice of the Julia set at any one time, how- 
ever, one can turn the view ninety degrees to view 
the two planes perpendicular to the current view- 
ing plane. One can also zoom infinitely into the 
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current image and roam vertically and horizontally 
in the current plane of view. Each HTML page is a 
virtual “chamber” of the Julia set with an address in 
(X,Y,Z,W) coordinates and a logical address in the 
oct-tree organization of the volume. 


When the user selects Approach the Oracle of 
(Chamber) ... the program invokes a Semantic 
Schema inspired by Cliff Pickover's book, Comput- 
ers and the Imagination: Visual Adventures 
Beyond the Edge, (New York: St. Martin's Press, 
1991: pp. 317-320). The dictionary (uncensored) 
contains references to Mathematics and Tantric 
practises. 


The Oracle invites the Explorer to “Inscribe (his/ 
her) Wisdom Upon the Chamber...” and instructs 
the Explorer to “Heed the Wisdom of Those who 
Have Passed this way”, which the Oracle shows to 
the Explorer. The Explorer is asked to “Identify” 
his/herself to the Oracle, to establish his/her per- 
sona in the MUD. 
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Fractal MUD is a catalyst for collaborative thought 
on the nature of Cyberspace — human social inter- 
action in a computer network — in which the 
computer is doing what computers do best: con- 
cretizing complex abstraction. 


The on-line documentation of Fractal MUD 
reminds the viewer that what the World-Wide Web 
can currently show is two dimensions removed 
from the actuality of the Quaternion Julia Set. 
There is a hot-link to Stewart Dickson's images of a 
Julia set pixel volume which has been reproduced 
in Polycarbonate using the Selective Laser Sintering 
technology at the Department of Mechanical Engi- 
neering, University of Texas at Austin. (See next 
article) 


We are thus encouraged to seek out physical arti- 
facts of Cyberspace which bring us closer to truth. 


Direct Sculpture of Deterministic 
Fractals 
Stewart Dickson 


On-Line Documentation of Direct Sculpture of 
Deterministic Fractals is available on the Internet 
(see reference [Ref3]). 


The Selective Laser Sintering (SLS) process has the 
ability to build on a particle-by-particle basis a 
computer-specified object as if it were a pixel vol- 
ume treated as a series of horizontally-oriented 
one-bit-per-pixel image slices. (See [Ref4] for a full 
description of the technology). The only require- 
ment to obtaining a solid final object is some 
abstract verification that the pixel volume is con- 
nected. Otherwise, there is no technical barrier to 
prevent building a disconnected or fractal volume. 
SLS is one of the only commercially available pro- 
cesses which can build this kind of object. 


The process was invented by Professor Joseph 
Beaman, Department of Mechanical Engineering, 
University of Texas at Austin. The commercial SLS 
machines are marketed by DTM Corporation, Aus- 
tin, Texas. They have a larger build envelope (a 
cylindrical diameter 12” by 14” axial) and build to 
higher precision than do the research machines at 
UT Austin, which build to roughly a 4” cube. 


However, the research machines, are open systems 
which routinely use custom software. The com- 
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mercial machines only accept a data interface in 
the form of the so-called ‘STL’ object file format, a 
polygon-based specification which is inappropriate 
for treating pixel volume data. UT Austin has 
defined a so-called “toggle” interface which is very 
nearly a one-bit-per-voxel run-length-encoding of 
laser scan-line specifications. I have calculated a 
nearly 25:1 data compression using the UT toggle 
file compared to the best that could be done with 
‘STL’. 

To date I have delivered nine data sets to UT Aus- 
tin for construction. The data sets are pixel vol- 
umes representing three cubes of decreasing 
parametric dimension (but constant pixel resolu- 
tion) of different iterative procedures: 


Z— Z + (-0.1103, —0.67037, 0, 0) 


Z- ez Z + e? (-0.745, 0.113, 0.01, 0.01) l 
and 


si ITC 
Z= e? Z + e ? (0.2809, —0.53,0,0) , where Z 
is a complex quaternion four-vector of the form 
(a, bi, cj, dk) . 


The pixel volume represents the three-space 
spanned by (1,470) (After V. A. Norton). 


Complete documentation of the project and pro- 
cess can be found on the Internet at World-Wide 
Web URL [Ref5]. 


The data sets I have created represent three-step 
“fractal zooms” in physical 3-D on three different 
procedures. The process can be extended up to at 
least an eight-inch cube at 100 to 200 voxels per 
inch, if the toggle-format interface ever becomes 
available on the larger, commercial SLS machines. 
Further, the same maximum resolution can be cal- 
culated on ever smaller regions of the domain over 
which the Julia set is defined, resulting in the reso- 
lution of ever finer detail. As you know, the pro- 
cess is bounded only by the floating-point 
numerical resolution of the machine on which the 
calculation is done. 


Selective Laser Sintering operates on theoretically 
any powdered thermoplastic material. Those cur- 
rently in use include polycarbonate, Nylon and 
machining wax. There is research currently under- 


way to build directly in metals. I propose to make 
fine art castings in Titanium of wax originals built 
using SLS. When the market becomes ripe for these 
objects, they will be available through The Will- 
iams Gallery ([Ref6)). 


The contact for the SLS portion of this project is 
Richard Crawford, Assistant Professor, Department 
of Mechanical Engineering, ETC 5.160, University 
of Texas at Austin, Austin, Texas 78712, (512)471- 
3030 (see [Ref7]). 


References 


e-mail: <dickson@mathart . com> 
Phone: (805) 494-6713 / FAX: (805) 371- 
1214 

URL: http://mathart.com or its mirror: 
http://www.wri.com/-mathart 
http://mathart.com/FractalMUD/ 
FractalMUD_home.html 

http: / /mathart.com/porttoli10/7 
SPD_Frac_portfolio.html or its T-1 
mirror Site, 

http: //www.wri.com/~mathart/port- 
folrto/SPD. Prac porttollo: html 
http: //ww.dtm-corp.com 
http://mathart .com and its mirror 
site http://www.wri.com/-mathart 
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lery/welcome.html 

e-mail: <wmgallery@aol.com> 

e-mail: <rhc@bian- 

chi .me.utexas.edu>. See also: http: / 
/sffoffice.me.utexas.edu/ 


[Ref7] 


Internet Postings 


Differential Fractal Geometry 
Dwayne wrote: 


I remember reading an article in the journal of clas- 
sical and quantum gravity a few years ago (late 
eighties?) by someone whose last name was some- 
thing like Nattle. I do remember that he was at a 
French University. He used assumed fractal dimen- 
sions for space-time co-ordinates, defined a fractal 
derivative, and claimed this automatically led to 
quantization of space-time. Sorry, I can't remember 
his exact name or what year it was. 


Laurent Nottale (this is his name) has presented a 
short review of his ideas at the FRACTAL 95 con- 
ference held in Marseille. He wrote a book two 
years ago about this problem (Fractal space-time 
and microphysics: towards a theory of scale relativ- 
ity, World Scientific, 1993). He told me that the first 
two chapters are absolutely easy to read (and are 
without any formulas). Try it! 


His e-mail is <nottale@mesiob.obspm. fr>. 


Locally Connected? 


In the sci.fractal FAQ, under “Q6h: Is the Mandel- 
brot set connected?” it says: 


Some connectedness definitions: 


Locally connected: X is locally connected if for 
every point p in X, for every open set U containing 
p, there is an open set V containing p and con- 
tained in the connected component of p in U. Le. 
X is locally connected if every connected compo- 
nent of every open subset is open in X. 


Could someone give examples of locally con- 
nected and not locally connected sets? (the defini- 
tion makes no intuitive sense to me). 

David Sharp 


Some Examples 


A connected but not locally connected set: 


Take the set which is the graph of sin(1/x) for pos- 
itive x. This will approach the y-axis arbitrarily 
closely. Now add the y-axis to the set. The result 
will be connected since you cannot place a bound- 
ary between the graph and the axis. It will not be 
locally connected since every sufficiently-small 
neighborhood of any point on the y-axis will not 
be connected. 


A locally-connected but not connected set: 


Two discrete points. It is obviously not connected 
but the neighborhood of each point is just that 
point and a point is connected. 

Joseph Hertzlinger 


New technique for 256-color M-set pic- 
tures 


When making 256-color Mandelbrot Set pictures, I 
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have long been bothered by the difficulty of 
approaching anything like an equal distribution of 
colors. The problem is in the use of dwell bound- 
aries; even the best choice of these does not pro- 
duce 255 regions, each containing the same 
number of pixels. 


Now, using the distance estimator algorithm 
(described in The Science of Fractal Images), I can 
produce such pictures. I think they are interesting 
enough to be considered by other subscribers to 
this list. The idea is simple: first calculate the (float- 
ing point) minimum distance for each pixel. Sort 
these in descending order (excluding those for 
which the distance is zero). Then assign a distance 
to each color-zone boundary so that each zone 
contains an equal number of pixels (distances). 
Now, for each pixel, determine its color- zone. The 
amount of computation is quite feasible, on a 
486DX. 


I have a GIF picture (SPIKE.GIF, 1024 X 768, 
about 165K) which I’m happy to send via ftp if 
someone has a home for it. The picture may be 
viewed by FRACTINT, although its M-Set parame- 
ters are not in the file. FRACTINT can be used to 
edit the colors to values other than my choice. 
Mike Freeman <mfreeman@unixg.ubc.ca> 


~Newton’s method & Julia sets 


I read somewhere that if you use Newton’s method 
to solve for the roots of a third-degree polynomial, 
and the points (in the complex plane) that do not 
converge to a root are plotted, then the result is a 
normal, second-degree Julia set. So I’ve been try- 
ing this for some time now, but I can’t seem to find 
equations with enough points that don’t converge 
to roots. How can I find equations that are hard to 
solve with Newton’s method? The book said some- 
thing about testing the ‘barycenter’... but what 


exactly is a ‘barycenter’?! Or is this Julia set/New- 
ton’s method thing a myth? Please, someone help! 
Christina Inloes 


Take any third or higher degree polynomial, say 


1. Consider a triangle with vertices at (X,, y,) , 
(X, >), and (X3, Y3) . The barycenter of the 
triangle is the “average” point of the three: 
(Ua +0 +X) 73, (Y, +> + V3) 73) - 

— Editor 


23-1. Scan the complex plane and apply New- 
ton's method to each point. Color the point one of 
three different colors depending on which of the 
three roots the point converges to (and every point 
you pick will converge to one of the three roots). 


When you are finished, you will have mapped the 
plane with three colors. The boundary between 
these three regions is the Julia set of the iterative 
equation given by Newton's method. 


Note that there are an infinite number of points in 
the Julia set, but your chances of hitting on a point 
that is a member of it is essentially nil. Think of it 
this way: you know that there is a real number 


equal to J2, but if you pick points from the real 
line at random and square them, you will never get 
a result of two. An even better example would be 


the circle whose radius is /2 (or two itself, for that 
matter). If you keep picking points at random you 
never hit one exactly on the circle. 


As you can get as close to 2 as you wish by suc- 
cessive approximation, you can get as close to the 
Julia set as you wish also. Just keep zooming in on 
the boundary. You find the points about the 
boundary take longer and longer to converge the 
closer you get. Mind you, they still converge but 
you can get arbitrarily large convergence periods 
by getting sufficiently close to the Julia set. 


Hope I didn’t misunderstand your question, and 
that this clears it up. 
Dennis C. De Mars 


HOP: Fractals in Motion 


HOP is a new fractal-generating program by 
Michael Peters and Randy Scott. HOP version 1.2 is 
now available. 


HOP is different that other fractal programs you 
may be familiar with. It generates “hopalong” frac- 
tals, so named because of their dynamic nature. 
Hopalong fractals differ from the better-known 
Mandelbrot fractals in a number of ways. The Man- 
delbrot set is based on complex numbers; hopal- 
ong fractals are based on real numbers. As 
opposed to the Mandelbrot set, which produces 
Static images, Hopalong fractals are moving 
images. Hopalong images are not infinite in scope 
(i.e., you cannot “zoom in” infinitely as you can on 
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a Mandelbrot image); rather, they are infinite in 
motion. Whereas the Mandelbrot set often pro- 
duces patterns that repeat in scope (as you “zoom 
in”) with slight variation, HOP often produces 
images that repeat, with slight variation, over time 
as the fractal develops dynamically. 


HOP Features 

e Fractint-like parameter files 

e Imports/exports 

e Fractint-compatible MAP palette files 

e MAP palette editor (registered version only) 


e GIF read/write screensaver for DOS, Windows, 
and OS/2 


e Resolutions up to 1280 x 1024 x 256 colors 
(unregistered version limited to 640 x 480 x 256 
colors) 


HOP Requirements 


IBM PC compatible computer w/640k RAM 
e Color VGA monitor 

e SVGA recommended for registered version 
e 386 or faster processor recommended 

e Math coprocessor recommended 


HOP is available at 
e FTP:rever.nmsu.edu (/pub/hop directory) 


e World Wide Web: http://rever.nmsu.edu/-ras/ 
hop 


e BBS:Fidonet 1:305/101 (505-646-2868; 14.4kbps 
supported) 


e CompuServe:GRAPHDEV forum (section 4: Frac- 
tal Sources) 
Randy A. Scott 


Fractal exploration (1) 


This may sound a little strange, but does anybody 
actually talk about fractals in this group? I'm a new- 
comer, and so far I’ve seen an educated debate 
about the appropriateness of posting binaries here, 
a lot of small talk about landscapes (try rec.hik- 
ing), and some interesting talk about video feed- 
back. Also requests for 24 bit fractal generators. 
What?! 


Does this mean that people out there have 


exhausted programs like Fractint? You’re done, 
explored it all, bored? Figured it out, blase? 
Hmmm. I know nothing about 24 bit fractals, and 
I'm sure they're really cool, but I don’t even want 
to get into it. It seems to have been established that 
you can color cycle the things, so why even 
bother? More is not necessarily better, and it seems 
ridiculous to want more of something when we 
have yet to master what we do have. 


Obviously the Mandelbrot set is bigger than any of 
us can hope to conquer. But we do make inroads, 
and we grasp certain underlying principles of orga- 
nization that appear all over, and at different 
scales. 


This for me is the beauty of fractalism, as much as 
the visual part. I think the best way to attack a 
problem like this is cooperatively, and this sort of a 
newsgroup would be ideally suited as a forum for 
sharing discoveries. Surely there are some fractal 
explorers out there. There must be! 


Has anybody ever attempted to make some kind of 
general vocabulary of the structures of the M set? 
Ever since they first spoke of the Seahorse Valley 
—and probably before — people have been see- 
ing organic creatures in there. I have my own jar- 
gon for these things and P'd be curious to know 
how others refer to these objects and areas. (Free 
color maps to the person who knows where to 
find the elephant-fungus (elephungus for short), or 
who knows how to cross the Butt-Brain Barrier!) 


There are huge ranges of topics for discussion, but 
a vocabulary seems like a logical place to start. 
Any interest? Anyone? Am I all alone here with my 
spirals? 

Jonathan N. Wolfe 


Fractal exploration (2) 


Ok, so somebody wants to discuss fractals! Very 
good. I too am trying to develop an understanding 
of the periodicities that occur in the M-set. 


Adrien Douady has a good article in The Beauty of 
Fractals that is hard, but worth going through. I 
can more or less understand the pattern of the 
orbits, and the “O” command in Fractint is really 
the ultimate fractal learning tool. What I have a 
harder time grasping is why that same periodicity 
that you see in the orbits manifests itself in the 
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Julia set of that point. 


It sorta makes sense, but is also miraculous. It 
allows you to go find a Julia set with an arbitrary 
number of arms, just by examining the orbit struc- 
ture along the inside of the main body of the M- 
set, and the degree of twistedness that you find in 
the orbits corresponds to lateral distance from the 
bulb on the set that is “influencing” the orbits. 


E.g. if you’re right under the period three bulb, the 
orbit shows three radial arms, and starts twisting as 
you go either right or left. The same thing happens 
in the Julia set, and most wonderfully of all, it hap- 
pens in the fine structure of M. 


Does anybody have an intuitive grasp of why this 
is? The structural relationship between the fine 
detail of the Mset and the local Julia sets seems to 
me to be the deepest, most stunning thing of all. 
One is a map in the C-plane, and one in the x- 
plane. Why should they mimic each other? How 
could this possibly be!? 


This is one of those pointless, unanswerable ques- 
tions, but if anybody has any thoughts, rumina- 
tions, meditations, or insights of any kind, do share 
them. 


Another possible topic to discuss: Has anybody 
played with the real or imaginary perturbations of 
the M set? This is when you start iterating the equa- 


tion z— z + c from some point z other than 0. It 


makes for very weird fractals indeed... 
Jonathan N. Wolfe 


Hypercomplex Fractals 
Jonathan N. Wolfe 


Has anybody explored in any great depth the 
Hypercomplex Mandelbrot set? I’ve recently gotten 
sucked in, and am suffering from that psychologi- 
cal condition characterized by sleeplessness, dizzi- 
ness, manic exuberance, and despair over the 
unconquerable, multidimensional complexity. 
(This may be a common state in this group, I don’t 
know.) 


I am starting to understand some of the patterns of 
organization of this 4D set, and am wondering if 
other people have insights into how it works too. 


All Pve been exploring is the “Sqr” function of the 
“Hypercomplex” fractal in Fractint and associated 


Hypercomplex Julia sets.” 


Very slight displacements of the j or k parameter 
lead to wonderful, ghostly contours, as the fractal 
fades in and out of the visible x-y plane. At least 
this is what I believe is going on. The coarse struc- 
ture of the set undergoes radical transformation 
with very small displacements in 7 or 7, but the real 
beauty is to be found down deep. There are, of 
course, plenty of “baby 'brots” embedded at all 
scales, and they too are 4-dimensional. Wait till you 
see a section of a 4D Seahorse Valley spiral vanish- 
ing and reappearing over and over as it spirals 
through the visible plane! 


I don’t think we'll be able really to grasp the multi- 
dimensionality until we can watch videos over and 
over of features of the set as ¡and Rare slowly var- 
ied. Pm working on this. 


Anybody have any thoughts on this set? 


Quartet Formula Development 
Noel Giffin 


I’ve had a few inquiries about some of the images I 
have posted to the alt. binaries.pic- 
tures.fractals newsgroup lately Some of this 
has been posted here before, but I have made 
improvements to the formula making it more gen- 
eral and also improving the bailout conditions. | 
also think my understanding is a little clearer now 
so pardon me if I reiterate. That is what we do best 
in this group anyway. My focus lately has been on 
the quartet fractal and it is responsible for most of 
my recent images. Someone asked if if was related 
to the Mandelbrot. 

The quartet fractal isn’t all that close to the Man- 
delbrot although, I admit it looks that way from 
some of the images. It is similar in some ways. It is 
a complex deterministic formula using escape time 
criteria for bailout but that is where the similarity 
stops. It’s another in a series of fractal formulae 
that I call inductive fractals due to an inductive 
self-reference in the iterative loop. The quartet 
fractal is actually closer to the lambda-sin fractal 


2. I think the “Hypercomplex Julia sets” 
Jonathan is discussing are the same as the 
quaternion Julia sets that Stewart Dickson dis- 
cusses elsewhere in this issue. 

— Editor 
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shown in The Science of Fractal Images, around 
page 160-162. There is a definite relationship 
between them, although I haven't proved it mathe- 
matically. The lambda-sin formula gives a repeat- 
ing fractal outline while my quartet formula gives a 
similar outline but fills the inside with very interest- 
ing fractal detail. This is due to the inversion which 
constantly turns the dynamics upside-down like 
turning an hourglass. The effect on the fractal 
image is dramatic. It no longer has an inside or an 
outside but shows a fractal that approaches a limit- 
set which I believe is the border of the lambda-sin 
function. Varying the initial conditions of the 
lambda-sin causes the fractal outline to skew and 
this happens in formula as well. What is fascinating 
is watching the internal structure change with dif- 
ferent initial parameters. The set seems to be an 
infinite extension of a repeating fractal geometry 
along the real axis. I haven’t checked this very far 
as fractint only has a max real ordinate of +32 or 
so, but it seems logical that this would be true from 
looking within those limits. tr is a full period. 

It isn’t that rewarding looking beyond a single 
period as the set is repeating. If there are any vari- 
ations, they are very minor. It is more interesting to 
change the initial conditions and observe a new 
geometry. It is essentially choosing a new julia set 
and the variations can be great. The trick is finding 
starting parameters that give good fractal detail. 
The P1 parameter is the starting point and can take 
on any complex value but the most interesting 
ones tend to be less than 1.0 (both real and imagi- 
nary) Actually, some of the most interesting are as 
these values approach zero. 


There seems to be a critical set of numbers where 

values less than these will produce image sets with 
large areas hitting the max iteration count. Num- 
bers slightly larger produce very intricate fractals. 
Any ideas on the shape of the critical set? I think it 
may either be points from the border of the 
lambda-sin function or maybe the Mandelbrot 
although the lambda-sin seems more likely. I 
haven't got enough points to really be sure of any- 
thing here. This is all very speculative. 

The P2 parameter is used to control one of the 
bailout conditions. Because I use a scaling inver- 
sion factor (P1), in the equation, and I have a pair 
of complex variables to control, I have to keep 
tabs on both of them or one could go to infinity 


while the other gets really small. They tend to see- 
saw back and forth due to the P1 inverted scaling, 
but if this number is less than 1.0 then the period 
can be many iterations and this gives the other 
variable a good chance to overflow. When the 
equation gets chaotic these values change drasti- 
cally. That’s why one of the bailouts is a large 
number approximately 1/2 of my floating maxi- 
mum value. The other bailout controls the fractal 
detail. Normally I use 64.0, but any value greater 
than zero can be used. Note that fractint defaults 
this value to zero and this is exactly wrong for this 
formula. You can see that this parameter controls 
the size of the disks that appear in the fractal 
periphery. A higher value has better detail (smaller 
disks) but takes longer to compute. The max itera- 
tion bailout can very anywhere from 64 to 32768 or 
whatever your max integer can be and as you can 
imagine it has the greatest impact on compute 
time, but some initial conditions, (those close to 
the critical values) require very high iteration limits 
to resolve the fractal. 

Here is the formula preceded by a couple of 
related ones. 

The quartet formula is a variation on an earlier 
fractal equation of mine called the quatro fractal. 

(named from its four point geometry) 

quatro(XYAXIS) { 

; Forget the squared term and iterate a trig function 

z=pixel, zp = (0,0): 
temp = z 
z = sin(z) - zp 
zp = temp, |zpl <= 4} 

It combines that formula with a variant which 
introduces the scaled inversion. The first formula 
of this type was the HT formula. 

ht(XAXIS) = { 

; Ah! You can use a variable in the inversion but 
now check for overflow. 

; Good results are found when the real 

ZE LO 
zp = pl/temp, (Izpl <= 64) && (Iz1 <= huge) 
} 

I was so happy with this fractal that I tried the 
technique on many of my earlier formulas and 
turned up the quartet formula which is very fasci- 
nating indeed. The quartet formula is: 
quartet { 

; The quatro-HT variation. 
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- Adding the inversion causes what looks like an 
infinite extension along the real axis. These should 
give a reasonable indication of the variations that 
this formula is capable of. 


Fractal Drawing Styles 
— Jesse Jones 


Ten different ways to color fractals are discussed. 
Many of these methods are discussed in The Sci- 
ence of Fractal Images [11, but the book can be dif- 
ficult to follow. Hopefully the material below will 
be easier to digest. Eight of these methods have 
been implemented in a program called Mandella 
for the Macintosh. 


GENERATING FRACTALS 


Typically fractals are generated using complex 
numbers. Every complex number can be repre- 
sented by a point on the plane. To generate a frac- 
tal we iterate some formula, for 

example z + z? + c. For each iterate we'll get a 
new z value. These points are called orbits. 

Here z is a point in the plane and cis a complex 
constant. When this is iterated one of two things 
can happen: z can blow up to infinity or z can 
become periodic. If z never goes to infinity it is 
inside the fractal. 

To simplify things we use a bailout value to see if z 
is headed for infinity, and a number called maxlter 
or hiCount. If the number of iterations exceeds 
maxlter we assume the point is periodic and there- 
fore in the fractal. 

After iterating a point we can color it based on the 
number of iterations or on the final z. Most fractal 
programs use the number of iterations to assign 
colors to points. This is called the “Level Set” 
method. 


THE DRAW STYLES 


(The constants defined below are special count 
values. Unity is used for draw styles that return a 
real number (say x). These numbers will be in the 
range 0 to 1. The count value is simply Unity x.) 


Constants 
Unity = 10000; 


mMaxCotint. = 65520: VE. .(* haCount <= max= 
Counce. *) 

inSet. = 65521;\t(* caught by period 
checker *) 

inFil = 65522;1t(* distance estimator 
filament *) 

nocount = 65535; 1€(*% no count value för 


Point *) 


Variables 
Z: starts as a point in the complex plane 
and is then iterated. 


orbits: an array that holds the value of 
z for each iteration. 


counT: the number of iterations. 


hiCount: the maximum number of itera- 
tions. 


bailout: the number used to determine 1f 
z has “blown up”. 


1. Level Sets 


This method uses the count value to index into the 


palette: 
index := (NoColors*count)/hiCount + 
FirstColor. 


However this will produce poor results. There are 
two problems with this approach: 1) If the image 
has been zoomed the lowest count in the image 
could be quite large. This means that a possibly 
large portion of the palette is never used. 2) The 
higher count values cluster close to the set. This 
causes a large amount of the palette to be dedi- 
cated to a very small area next to the fractal. This 
problem will become worse if the hiCount is 
increased. 
To alleviate these problems we can use three cut- 
offs: high, middle, and low. The high cutoff will be 
used for the hiCount. If the count is between the 
low and middle cutoffs we use the count value to 
index into the palette. if the count is greater then 
the middle cutoff we use the last color in the pal- 
ette. If the count is less then the low cutoff use 
some background color. Putting it all together we 
have: 
IF count < low THEN 

index := BackColor; 
ELSIF count >= high THEN 

index := SetColor; 
ELSIF count > middle THEN 

index := LastColor; 
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ELSE range := LONG(middle - low); 

delta := LONG(count - low); 

index := (delta*NoColors) DIV range + 
FirstColor; 

END; 


This is much better but we have only partially 
solved the problem with the counts clustering 
close to the fractal. To bring out the most detail we 
need to somehow use the bulk of the palette for 
the lower count values. One possibility is to use a 
formula like this: 
newCount. <= hiCount > 

(count /hiCount)*(1/root) 

This formula will take a count from 0 to hicount 
and return a new count value in the same range if 
root > 0. If root = 1 then newCount = count. If 
root > 1 then the low counts account for the 
majority of newCounts. This is what we want since 
the low counts will then be assigned to most of the 
palette. If root < 1 most of the newCounts come 
from the higher count values. This can sometimes 
be interesting. To integrate this with cutoffs we can 
write: 

range:= LONG(middle - low); 


delta := LONG(count - low); 

new := range*(delta/range)*(1/root) ; 
index := (new*NoColors) DIV range + 
FirstColor; 


Another idea is to allow the user to somehow 

draw a graph to control the count to palette map- 
ping. An easy and powerful method is to use Bez- 
ier curves. The December 1986 Byte has a good 
article on Bezier curves. The rest of the drawing 
styles return real numbers (say x). To index into 
the palette the real number must be somehow 
scaled into a number from 0 to 1. Then the palette 
index is simply: 

index := (NoColors-1)*scale + FirstColor 
To do the scaling we need to find the largest x. 
Then scale := x/maxX. In an interactive program 
maxX can sometimes be hard to find. 

Many programs using the Level Set method store 
each count value in an array. This array can then 
be used to redraw the image (perhaps after a 
change in the root scaling). To integrate this with 
the other drawing styles we can multiply scale by 
some large number Unity. Then 0..Unity will be 
inside the palette and Unity+1 will be inside the 
set. 
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2. ColorSet 


Here we color the fractal and leave everything 
outside the fractal black. To get a color we use the 
final z value. This style works well with Julia type 
fractals and period checking. The code looks like 
this: 

IF count >= hiCount THEN 


mad o= TEE 2 EAS 

scale := mag - TRUNC(mag); (* ignore 
integer part *) 

count := Unity*scale; 

ELSE 

count” Y=: 0. (* use FirstColor *) 


END; 


3. Decomposition 


In this style we treat z as a point in the plane and 
use the angle z makes with the x axis to index into 
the palette. In the Science of Fractal Images [1] they 
use “Binary Decomposition”, i.e. they use only two 
colors. We will use all the colors. 

IF count >= hiCount THEN 


count = Unity+L; 
ELSE 
theta. w= AreTan(z.1/z-r) + PI/2; 
scale := theta/PI; 
count := Unity*scale; 
END; 
4. Stalks 


This method is due to Pickover. I found it in Com- 
puters and the Imagination page 170. The idea is 
to color a point differently if either the real part or 
imaginary part of an orbit are less then a small 
number (e.g. 0.01). The result is a standard looking 
image but with tendrils of varying widths and col- 
ors going every which way. The code I use in Man- 
della looks like this: 

VAR 

minz REAL; (* ABS value of smallest 

component in each orbit *) 

epsilon: REAL; (* A small number entered 

by the: user, *) IF count < hiCount “THEN 
IF minZ < epsilon THEN 


scale := minZ/epsilon 
count := loCount + scale* (midCount - 
loCount); 


END; 

END; 
Note that the interior of the set could also be col- 
ored so that it has stalks. 


5. Continuous Potential 


Here we use the magnitude of z and the count 
value to get a smoother image. The potential is cal- 
culated like this: 

Pot. += LR(Mad(z2)1)/ (2 count.) 
To scale the potential into the range [0..1] we 
need to find the largest potential: | 


maxPot := Ln(bailout)/(2*1loCount) 
Then 
scale := pot/maxPot 


As count increases scale will go rapidly to zero. To 
smooth things out we will use a variable called 
slope. Using a large bailout value will also help 
matters (500 works well for the classic Mandelbrot 
fractal). The code follows: 

IF count s= -hiCount THEN 


count. s= Unicy+l> 

ELSE 

pot. += EniMagtz))t2 Count): 

maxPot := Ln(bailout)/(2*low); 

scale := pot/maxPot; 

scale := scale” (1/slope); 

scale := 1 - scale; (* shift high counts 
to the right of palette *) 

count := Unity*scale; 

END; 


6. Distance Estimator 


The distance estimator method is used to identify 
points that are close to the set. The user inputs an 
epsilon that determines how close to the set points 
have to be before they are drawn (0.6 works fairly 
well for the entire Mandelbrot Set). As the user 
zooms in epsilon should be decreased. 

Pages 196-198 of The Science of Fractal Images [1] 
discuss the mathematics behind the distance esti- 
mator method. There are also some striking black 
and white illustrations. I wasn’t able to get any- 
thing that looked as good. I suspect the higher res- 
olutions used in the book makes a big difference 
with this method. 

In Mandella 6.0 the Distance Estimator draws 
points close to the set using a special color. Points 


farther away are drawn using the LevelSet method. 
Earlier versions used the distance to color points 
far from the set but I was never happy with the 
results. 

overflow ts 1L0E32% 

delta := epsilon*dx; (* where dx is dis- 

tance between pixels *) 

IF count. >= hiCount “THEN 


count. ie UNTIL 
ELSE 

¿Primesr x= 0; 
zPrimest <= 0; 

Kk ss Ty 


WHILE (k <= count) AND (ABS(zPrime.r) < 
overflow) AND (ABS(zPrime.i) < over- 
flow) DO 


zPrime := Mult(orbits[k], 
2Prime.r se. 2*zPrime.r + i; 


zPrime) ; 


zPrime.i := 2*zPrime.i; 


INC (k) 

END; 

IF k <= count THEN 

count = inFil; (~ point is very close 
to set *) 

ELSE 

magl <= Sgqrt(Mág(tz))5 

mag2 := Sqrt (Mag(zPrime)); 

dist := Ln(magl)*(magl/mag2); 

IF dist < delta THEN 

count: += Inds. (* point 1s fairly close 
to set *) 

ELSE 


(* point is far from set: treat like 
LevelSet *) 


END; 
END; 
END; 


7. Biomorphs 


I found this method on page 104 of Computers, 
Pattern, Chaos and Beauty [3]. It's a simple method 
that works well with low iteration levels. The fig- 
ures produced often resemble invertebrate life- 
forms. It works by setting the color of a point to a 
special value if the real and imaginary coordinate 
of z are greater then the bailout value. Otherwise 
the point is colored using the LevelSet method. 

IF count < hiCount THEN 
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IF (ABS(z.r) < bailout) OR (ABS(z.1) < 
bailout) THEN 

(* treat like LevelSet *) 

ELSE 

count := inFil; (* use the filament 
color *) 

END; 

END; 


8) Smurf 


This method is due to the Brainy Smurf (aka Rich- 
ard Hughes). The idea is to use either the real or 
imaginary components of z to color the image. The 
code looks like this: 

IF count < hiCount THEN 


gap.r := Size.r/FLOAT(portRect.right - 
portRect.left) ; 

gap.i := Size.i/FLOAT(portRect.bottom - 
portRect.top) ; 

toler := MinReal(gap.r, gap.1)/2.0; 
END; 


END InitPeriodChecker; PROCEDURE Equal 
(x, y: REAL): BOOLEAN; 

BEGIN 

RETURN ABS(x - y) < toler; 

END Equal; 


PROCEDURE DoPeriod (z: Complex); (* Call this 
for each orbit *) 


BEGIN 

IF WasInSet THEN 

IF Equal(z.r, lastZ.r) AND Equal(z.i, lastZ.i) THEN 
oldCount := count; 


scale := bailout/ABS(z.r); (* get a num- 
ber from 0 to 1 *) | count := inSet; 
scale := scale” (1/4); (* force the num- ELSE 
ber closer to 1 *) INC(Step); 
count := count*scale; (* adjust count *) IF Step > maxStep THEN 
END; step <=. 1; 
maxStep := 2*maxStep; 
lastZ := Z; 
9. Periodic END: 
This one works only if you are using period END; 
checking. The idea is to color the set based on the END; 
number of iterations it takes before a point END DoPeriod; The Periodic draw style code is 
becomes periodic. The period checker I use is then: 
listed below. It’s designed to quickly catch the IF count = inSet THEN 
points that rapidly become periodic. count := oldCount; 
VAR ELSIF count >= hiCount THEN 
WasInSet: BOOLEAN; (* TRUE if the last count := hiCount - 1; (* Use LastColor *) 
point was in the set *) ELSE 
toler : REAL; (* x = y if ABS(x - y) < count := 0; (* Use FirstColor *) 
toler *) END: 
Step : INTEGER; (* current period length 
e) 
maxStep : INTEGER; (* maximum period to 10. Ratios 
check for *) This is a simple method that turns out to work 
lastz : Complex; (* z value to compare rather well. The results are similar to decomposi- 
with *) tion but usually more interesting. The idea is to use 
oldCount: CARDINAL; (* unsigned integer the ratio of z.r and z.i to color the image. The code 
A a 
IF count < hiCount THEN 
VAT gaDt ROME max := Max(z.r, z.i); 
PEGEEN min := Min(z.r, z.i); 
IF WasInSet THEN scale := min/max; 
ci Sri count := Unity*scale; 
Step := maxStep; lastZ.r := 0.0; END: 
Tastístk t= 0203 
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ADDRESS 


Jesse Jones 
Usenet: jesjones@milton.u.washington.edu 
CServe: 73627,152 


Submitting Articles 
Here are some guidelines for submitting articles for 
publication in Amygdala. 


1. Type your article as you would like to see it 
appear in the newsletter. Please do not send hand- 
written drafts! 


Send it (my order of preference): (A) On 3.5” dis- 
kette in Macintosh format for FrameMaker, 
MacWrite, WriteNow, or MS Word; or as a text file. 
Please also enclose paper copy so I can see your 


Amygdala is published by Rollo Silver, ten issues 
per year (hopefully). 


Address articles, letters, comments, subscription 
orders, to: Amygdala; Box 219; San Cristobal, 
NM 87564; USA. 

Email: rsilver@artvark.com 

Phone/FAX: 505/586-0197 

Compuserve: 71174,1453 

This issue was produced using FrameMaker on a 
Macintosh IIcx computer. 


intent. (B) On 3.5” diskette in IBM format, text file. 
I have no way to deal with 5.25” diskettes. (C) 
Paper copy. 


2. Illustrations should be either greyscale (suitable 
for halftoning) or black/white; not color! (A) Nor- 
mally, illustrations will be printed in full column 
width, so you should make them 3.25” wide, if 
possible — provided that they’re 300 dpi resolu- 
tion. If they’re grainier, make them larger if possi- 
ble, so that they’ll look good when reduced to 300 
dpi. (B) Make sure that you clearly indicate which 
illustrations go where in the text! (C) All in all, it’s 
better not to have captions welded into your pic- 
tures. Let me put them in ad lib. (D) I can handle 
illustrations on diskette in MacPaint, MacDraw, 
CricketDraw, Photoshop, or Adobe Illustrator for- 
mats. 


3. Please send along a short biographical note, 
which I will try to publish in the same issue as 
your article. 


4, Please include your telephone number, in case 
I have to reach you in a hurry with questions. 


Circulation/Renewal 


As of September 19, 1995, Amygdala has 391 sub- 
scribers. 


The address used for mailing this issue to you has 
a small notation at the lower right indicating when 
your subscription expires. 


A35 indicates that this is your last issue. A36 indi- 
cates that your subscription expires with the next 
issue. In either case, please use the form on the 
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Jewel Box 


Offset Litho Print 


We are offering an offset litho print of Rollo 
Silver's beautiful fractal image, Jewel Box, 
printed on 100 pound glossy archival stock. 
You can see a small replica of the image on 
the card to the left. 


With all the clarity and crispness the card 
possesses, you can imagine the stunning 
quality of a print ten times as large, with ten 
times the resolution! 


The print is 24 x 32", with the actual image 
area being 21.25 x 27": 8,232 x 6,480 pixels. 
To see the Jewel Box in a larger size and to 
see many more of Rollo Silver's fractal 
images, check out the web site: 


http: //www.artvark.com/artvark/ 


The regular price for the print is $29.95, but 
until October 15 we are offering it to 
Amygdala subscribers for $24.95. 


The print will be shipped to you in a stout mailing tube via Priority 
Mail (U.S.) or, outside the U.S, via Air Mail or Surface Mail at your 
choice. Bear in mind that Surface Mail can take many weeks to get 
to you. 

> COPIES OF JEWEL BOK, O $24,959 CA siii 

Postage & Handling (18 oz per item): 

© First Class/ Priority Mail in U.S. @ 93,00 Cd 

© Air Mail, Canada $3.50, Mexico $4.00, W. Hemisphere $6.00 .......... 

© Air Mail, Europe $8.00, Asia/Africa/Pacific Rim: $11.00.................... 

Q Surface Mail, Canada/Mexico $3.00, Others $5.00 woo... 


Payment enclosed (please indicate which): 
check d MasterCard Y VISA Q American Express J Other Q 


Credit Card +: Expiration Date: 
Name: 

Address: 

City, State, Postal Code, Country, etc.: 

Signature (for credit card charge): 


Telephone: 


Send order to AMYGDALA e Box 219 + San Cristobal, NM 87564 e USA 
Telephone: (505) 586-0197 
Email: rollo@artvark.com 


Please pay with a credit card, a check payable through a US bank, or money 
order, postal money order, bank draft, cash, etc. U.S. Funds, please! 


September 12, 1995 


JEWEL BOX 
Supergloss print 30" x 40" 
O Rollo Silver 1993 


RECURSIVE IMAGES 


Fractal Art by 


Rollo Silver 


Lama Round House e PO Box 219 
San Cristobal, NM 87564 
(505) 586-0197 


